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I. INTRODUCTION 

Completely-positive, trace-preserving (CPTP) quan¬ 
tum maps capture a vast diversity of quantum dynami¬ 
cal evolutions, including arbitrary open-system dynamics 
such as decoherence, measurement, and thermal relax¬ 
ation Consequently, the thermodynamic analysis 

of processes described by CPTP maps is a major issue 
in the development of quantum thermodynamics M- 
One of the main tools of such thermodynamic analysis 
are fluctuation theorems, since they reveal the statistical 
properties of thermodynamic quantities such as work or 
entrop 5 ^roduction along arbitrary nonequilibrium pro¬ 
cesses (8l-[l(]l|. Therefore, fluctuation theorems for arbi¬ 
trary CPTP maps will be relevant to understand the role 
of quantum effects in thermodynamics. 

In recent years, there have been several derivations of 
fluctuations theorems for specific classes of CPTP maps 
falling into two broad categories: detailed fluctuation 
theorems for quantum trajectories and fluctuation the¬ 
orems for thermodynamic variables, such as work and 
entropy. Campisi et al. obtained a detailed fluctuation 
theorem for a unitary, driven evolution punctuated by 
unital maps - maps for which the identity matrix is 
invariant, such as projective measurements - [H d. 
This work was followed up and extended by Watanabe et 
al. [l3l| . General quantum Markov semigroups were ex¬ 
plored by Crooks using time-reversed or dual maps 
which was then applied by Horowitz et al. to nonequi¬ 
librium quantum jump trajectories [id- An alter¬ 
native, operator formulation for driven Lindbald master 
equations was independently developed by Chetrite and 
Mallick [l^. Fluctuation theorems under unital CPTP 
maps for thermodynamic quantities, like work, energy 
and information-theoretic entropy, have appeared in nu¬ 
merous works did: while predictions for non-unital 
CPTP maps usually take the form of an integral fluctu¬ 
ation theorem with a so-called correction [isl - d • 

In this paper we present a general formalism based on 
a generalized detailed balance condition that includes ex¬ 


tends many of the previous results without the need for a 
correction term. As a consequence, our result clarifies the 
minimal hypotheses needed to derive a fluctuation theo¬ 
rem for quantum maps. Our theorem is independent of 
the physical nature of the process that induces the CPTP 
map. This is a relevant feature. It makes the fluctuation 
theorem general enough to be applied to situations far 
from equilibrium, like systems in contact with coherent 
reservoirs [l^ d- Moreover, such a general result 
could be useful to analyze the thermodynamics of quan¬ 
tum processes whose physical details are not completely 
known, such as decoherence or quantum collapse. 

The paper is organized as follows. In Sec.|n]we review 
the theory of CPTP maps and the Kraus representation, 
introducing the dual map, necessary to state the fluctu¬ 
ation theorem. In Sec. uni we prove the general theorem 
for single maps and for a series of concatenated maps. 
Some applications are discussed in Sec. |TVl Finally, in 
Sec.lv] we summarize our results and present the main 
conclusions of the paper. 


II. QUANTUM OPERATIONS AND DUAL 
DYNAMICS 

Consider a generic CPTP quantum map p ^ p' = £(p) 
acting on the density matrix p of a quantum system. Any 
CPTP map admits a Kraus representation in terms of a 
collection of linear operators {Mk} as Hi 

£:(p)=^ffc(p)=^MfcpA4^ (1) 

k k 

with = 1, a condition that ensures the trace¬ 

preserving property of the map £. It is important to 
stress that the choice of {Mk} is not unique: any family 
of operators M/ = J2k UikMk, with Uik the entries of a 
unitary matrix, is also a valid Kraus representation. Not 
even the number of Kraus operators is unique for a given 
map. For instance, if the Hilbert space of the system has 
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finite dimension N, there exists a Kraus representation 
for any map with at most iV^ operators. However, using 
more than N'^ operators is sometimes necessary for a 
complete description of the physical process associated 
to the map (as we will see below). 


A. Quantum trajectories and nonselective states 


The Kraus representation © is not just a mathemat¬ 
ical way of writing the map; it also provides a physical 
picture of the map as a random transformation of pure 
states. A specific representation decomposes the map 
into a number of operations Ski') = Each op¬ 

eration transforms a pure state into a new pure state 


Wk) = 


Mk\ii^) 

\\Mk\^)\ 


( 2 ) 


with probability pfc( IV')) = \\Mk |V') IP (EfePfc(IV')) = !)■ 
This picture extends to mixed states of the form p = 
'lliPi \'4’i) which represents a classical ensemble of 
pure states \’4>i) each sampled with probability pi. Thus, 
the probability that operation k occurs is 


Pkip) = = Tr[£:fc(p)] 

i 


( 3 ) 


and the final state conditioned on this operation is 

, _ Mk IV'i) (V'il ^k _ ^kip) 

Pk-2^P^ ||M,|v^,)|p -pkip)- 


( 4 ) 


If we know which operation Sk has occurred, then k can 
be seen as the outcome of a generalized measurement 
and p'p. as the selective post-measurement state of the sys¬ 
tem. If we do not know which operation took place (or 
we decide not to incorporate that information into our 
description), then the state after the transformation is 
p' = Sip) = J2kPkip)p'k^ usually referred to as the non¬ 
selective post-measurement state (although the transfor¬ 
mation given by the map p' = Sip) does not necessarily 
imply any measurement and not even a specific Kraus 
representation). This setup defines an efficient general¬ 
ized measurement in quantum mechanics, more restric¬ 
tive than generalized measurements where the observer 
has access only to a function fjk) of the operation index 
k, which may not be one-to-one |24| . 

A generic quantum evolution is described by a concate- 
nation of maps Sr with Kraus operators Mj, . For the 
initial state p(0), the nonselective state evolves as 


pir) = SrSr-l . ■ .SipiO). (5) 


This density matrix p(r) can be interpreted as the av¬ 
erage of the stochastic evolution. If the initial state is 
pure p(0) = |^(0))('!/'(0)|, then a stochastic trajectory 
7 = (fci, k 2 , ■ ■ - kr) is given by the operations that oc¬ 
curred in the application of map Sr and determines the 
evolution of the pure state: 

|^(r)) = |V>(0)). (6) 


B. Dual dynamics 


Now consider a particular Kraus representation of 
a map S = ^j^Sk-, and suppose that the map has 
a positive-definite invariant state tt (not necessarily 
unique), i.e.. Sin) = tt. For such maps, we introduce 
an auxiliary or dual map S with respect to tt and to a 
fixed, arbitrary unitary or anti-unitary operator A. In¬ 
spired by Crooks, we define this dual map through the 
equality 0[il 


Tr (tt)] 


Tr 


Skffik^i^) 


( 7 ) 


where tt = AnA'^ is the invariant state transformed by A. 
Equation © states that the probability of observing the 
outcome fci followed by ^2 when we apply the map twice 
to the invariant state tt equals the probability of observ¬ 
ing the reverse outcome —^2 followed by fci— when the 
dual map is applied twice to n. In this way, the dual 
map induces a dynamics in the invariant state that is the 
reverse of the original one. Following the derivation in¬ 
troduced by Crooks in Ref. [3, one can prove that the 
Kraus operators of the dual map are given by 

Mk = An^ A^. ( 8 ) 

Trace preservation ij^k ^l^k = 1) follows immediately 
from Sin) = tt, and one can verify that the dual map 
preserves the dual invariant state. Sin) = n. 

The inclusion of the operator A in the definition of the 
dual map is not mathematically necessary to derive the 
fluctuation theorem. In fact, A does not appear in the 
original definition by Crooks [l^ . However, in some situ¬ 
ations an appropriate choice of the operator A is needed 
to find a dual dynamics with a precise physical interpre¬ 
tation or that is suitable of being implemented in the 
laboratory. The customary choice is the time-reversal 
operator A = Q that changes the sign of odd variables, 
like linear and angular momenta. 0 is an anti-linear, 
anti-unitary operator, satisfying 0^ = 0^'0 = 00^' = 
1 im. For instance, 0 acts on a spinless particle by 
complex conjugation of the wave function in the position 
representation. The need of 0 in the definition of the 
dual process is clear, for example, if the map is a unitary 
evolution, i.e., a map given by a unique Kraus operator 
U with C/^ = U~^. In that case the invariant state is pro¬ 
portional to the identity matrix and the dual dynamics 
reads 


L( = 0[/t0t. (9) 

The dual map is again a unitary evolution given by the 
unitary operator U and corresponds to the operational 
time reversal of the original unitary evolution given by 
U [ 2 ^. For instance, if U is the evolution of a system 
under a constant Hamiltonian H,U = and H is 

time-reversal invariant, \H, 0] = 0, then JJ = U, i.e., the 
dual map is identical to the original one. On the other 
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hand, if the Hamiltonian depends on time according to 
some protocol, and U is the evolution between t = 0 
and t = T, then U is the evolution that results when the 
protocol is reversed (which is, in general, different from 
C/t). 

The operator A can also account for other transfor¬ 
mations of the system state that are necessary to exploit 
dynamical and static symmetries. In fact, this freedom 
has a classical counterpart in fluctuation theorems that 
incorporate various symmetry transformations [27l - l^ . 

III. FLUCTUATION THEOREMS 
A. Nonequilibrium potential and detailed balance 

We now prove a general fluctuation theorem for a large 
family of CPTP maps. To begin our introduction of these 
maps, let us focus on an important class of maps that 
admit the following Kraus representation 

Mji = Oji ItTj) (tt.I , (10) 

in terms of the eigenstates {|7ri)} of the invariant den¬ 
sity TT. Here the Kraus operators are labeled by two 
indices (i,j) that identify jumps or transitions between 
eigenstates of tt, \'Ki) —>■ \TTj), occurring with probability 
W^ji Ki) IP = These maps are special in that a 

single application of £ destroys any coherences between 
eigenstates of tt in the initial state p, reducing the sub¬ 
sequent action of the map to a classical Markov chain on 
the eigenstates {|7ri)}. Therefore, the dynamics induced 
by CPTP maps of the form (fTUl) is essentially classical. 
On the other hand, quantum effects arise if the Kraus 
operators are linear combinations of the transition oper¬ 
ators \TTj) (TTil, preserving coherences between eigenstates 
of the invariant density matrix. 

The family of maps that obey a fluctuation theorem go 
slightly beyond the “classical” case outlined above (nnj. 
To make this family precise, we assign to each eigenstate 
iTTi), whose strictly positive eigenvalue is denoted by 7r(i), 
a nonequillibrium potential, similar to the one used in the 
classical Hatano-Sasa theorem [s^, 

d>^(i) =-ln7r(i). (11) 

Then the maps that obey our fluctuation theorem are 
those where each Kraus operator Mk is formed from a 
superposition of jump operators, all of them inducing 
the same change in nonequilibrium potential A^T^{k): 

Mk=J2m'^i\TTj){TTi\, (12) 

with = 0 if <l’ 7 r(j) — 'h 7 i-(*) ^ A$^(fc). That is, 
by measuring the operation Mk we know without uncer¬ 
tainty the change in the nonequilibrium potential, even 
though that change could have occurred through a su¬ 
perposition of jumps. One example of this construction 


is a harmonic oscillator coupled to an equilibrium reser¬ 
voir of resonant photons at temperature T [sij . Here, the 
nonequilibrium potential is the energy of each eigenstate, 
divided by kT, and the change in the nonequilibrium po¬ 
tential in a transition is proportional to the energy trans¬ 
ferred to the reservoir of photons as heat. By measuring 
the reservoir we are able to detect jumps in the oscilla¬ 
tor, but the measurement, in general, does not provide 
information about the system state |3l| . 

It is straightforward to check that condition (USD is 
equivalent to 

[Mfe,ln7r] = A$^(fc)Mfc 

[Mllmr] = -A^^{k)Ml (13) 

and, consequently [M^Mfc,ln7r] = = 0. These 

commutation relations are similar to those satisfied by 
the Lindblad operators that appear in Davies’ theory of 
systems weakly coupled to thermal baths (see below and 
They indicate that the pair Mk, M^ acts as lad¬ 
der operators, inducing jumps between the eigenstates 
\'Ki) of TT with a fixed change A$,n.(fc) in the nonequi¬ 
librium potential $. Finally, m ensures that the dual 
Kraus operators obey a generalized detailed balance con¬ 
dition 

Mk = AmIa'’ (14) 

that can be obtained by plugging m into ([5|). One can 
also prove that the form (IT^ is the only one for which 
the dual operators Mk in ([5|) are proportional to AM^A^. 
Remarkably, for maps with multiple invariant states the 
A<i>^(fc) do not depend on the specific invariant state tt 
chosen to define the nonequilibrium potential and the 
dual dynamics [s^. In other words, the set of values 
A4>,n.(fc) is a property of the map £. 

B. Fluctuation theorem for a single CPTP map 

The basis of our fluctuation theorem is codified in the 
proportionally between Kraus operators and their dual 
counterpart in (ED- This generalized detailed balance 
condition connects the probability to observe a given 
jump, say k, with the probability to observe the same 
jump in the dual dynamics. Specifically, suppose that we 
initially prepare the system in the pure state |V’n), and 
then apply the map £, registering the occurrence of the 
operation k. We then perform a quantum yes/no mea¬ 
surement of the projector \(j)m) The subscripts n 

and m are added to the initial and final states so that 
later on we can consider measurements of arbitrary ob¬ 
servables with eigenstates |^„) and \4>rn)- 

Now, let p{m, k\n) be the probability that given an 
initial state \ipn) we observe operation k and the final 
state that is, the probability to observe the jump 

IV'n) —t \(t>m) under the action of Mk- Let p{n,k\m) be 
the probability to observe the inverse jump \(j)m) —t \4’n), 
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with [tp) = A\ip)^ under the action of the dual opera¬ 
tor Mi~. Using (HI, the ratio of these two conditional 
probabilities is 


p{m, k\n) _ I {ipn) p _ | {4>m\Mk li’n) P 

p{n,k\m) \{^n\Mk\pm)\'^ \{lpn\APMkA\pm)\'^ 

_ I {4’m \ Mk \'4>n) I _ 1 _ -A$,(fe) /-i r\ 

" I !</,„) |2 - 

Equation (ITSl) can be considered as a modified detailed 
balance relation for the operation £k and its dual £k, 
which remarkably is independent of the initial and final 
states. 

Suppose now that we prepare the system in the initial 
mixture pi = '^^Piiji) \ipn) (t/'nl and apply the map £. 
By measuring the initial state |'0n), the operation £k and 
a final state \4>m) we obtain a trajectory (m, k, n) that is 
observed with a probability pirn, k, n) = p{m, k\n)pi{n). 
We compare this to a dual process induced by the map 
£ applied to the initial state pi = \4>m) {4>m\- 

The dual trajectory (n, k,m) is given as well by the ini¬ 
tial state \4>m), the dual operation £k and the final state 
\ipn), and it is observed with probability p{n,k,m) = 
p{n,k\m)pf{m). The ratio of the probability to observe 
a trajectory (n, k, m) and the probability to observe the 
reverse trajectory (m, fc, n) in the dual process is then, 
from (fTSl) . 


thermal baths, etc. A trajectory now is given by the ini¬ 
tial |'0n) and final states \4>m) and the outcomes kr of all 
the measurements associated to the maps r = 1,2,... ,R: 
7 = (n, fci, ^ 2 ,..., fc_R, to). Each map £r has a Kraus rep- 
resentation, given by the operators , and an invari¬ 
ant state 7r(’’^ for which the dual map £r and the nonequi¬ 
librium potential $^(r) (i) are defined as in Eqs. dS]) and 

(HU. 

To derive the fluctuation theorem, we reverse the con¬ 
catenation of maps. We define the dual process as 
Cl = £i... £r ... £r-i£r (notice that, for i? > 1, in gen¬ 
eral, Cl ^Cl, i.e., the dual process does not coincide with 
the dual map of U). If each map obeys condition (fT^ [or, 
equivalently, (HI], then we get the following symmetry 
relation 


p{m,kR,... ,ki\n) 
p{n,ki,... ,kR\m) 




exp 


R 

(kr) 

r—1 


(19) 


A detailed fluctuation theorem can be now obtained 
by comparing the probability of a trajectory 7 = 
(n, ki,..., kR, to) in the forward process and the prob¬ 
ability of the inverse trajectory 7 = (to, fe/j ,... ,ki,n) in 
the dual process: 


S(n, k, to) = In = a(n, to) — A4>^(fc), (16) 

p(m, k, n) 

where a(n,'m) = — Inpf (to)-|-I npi(n) is a boundary 
term, only depending on the initial state of the process 
Pi and the initial state of the dual pi. The quantity E 
is a measure of how different the original and the dual 
trajectories are. In particular, when the dual is the time 
reversed process (see below), E is a measure of the irre¬ 
versibility of the process for a given trajectory. In the rest 
of the paper we will show that it can be identified with 
an entropy production in many situations of interest. 

A Jarzynski-type intergral fluctuation theorem imme¬ 
diately follows from (fT^ : 

(e-^) = l, (17) 

where the average is over forward trajectories, p(n, k, to). 
Finally by Jensen’s inequality (e^) > we have the 
second-law-like inequality 

(E) = (a) - (A$,) > 0. (18) 


C. Fluctuation theorem for concatenated maps 

Our fluctuation theorems dmHzi) can be easily ex¬ 
tended to a concatenation of CPTP maps, £l = 
£r£r-i ... £r ■. ■ £1, which is the case of general Markov 
quantum evolution, unitary evolution punctuated by pro¬ 
jective measurements, driven systems in contact with 


E( 7 ) = In 


Pjl) 

P{l) 


R 

a(n,m) - ^ A$^M(fcr), 

r—1 


( 20 ) 


with a corresponding integral fluctuation theorem that 
follows readily, like in (ITTl) . Thus, for a concatenation of 
maps implemented in sequence, we merely have to add 
the changes in the nonequilibrium potential along the 
trajectory. Notice also that we effectively used a Kraus 
representation for the map £1 where each Kraus opera¬ 
tor was labeled with the sequence {fci,..., fe/j}, requiring 
possibly many more than the necessary operators. 

A clear interpretation of E( 7 ) arises if we consider the 
concatenation of the same map £, acting on the station¬ 
ary density matrix tt, and the corresponding dual process 
acting on ff. In this case Pi(n) = tt ^ and Pf(m) = tt ^, 
yielding 

R 

E(7) = ln7r„ - InTTm - ^ A$^(A:r) = 0 (21) 

r—1 


for any trajectory 7 . This is expected from the (modified) 
Crooks dehnition 0 : the original and the dual maps act¬ 
ing on TT and if, respectively, produce a trajectory 7 and 
its reverse 7 with identical probability. Therefore, E can 
be considered as a measure of the distinguishability of 
the original and the dual process, but also as a measure 
of how far the system is from the stationary state. These 
two equivalent interpretations are familiar in thermody¬ 
namics when TT is an equilibrium state: the dual is the 
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reverse process and £( 7 ) is the entropy production which 
measures both irreversibility and departure from equilib¬ 
rium In more general situations, S(j) is the part of 
the entropy production due to the fact that the state of 
the system does not coincide with the stationary state. 
This can occur in the transient from a nonsteady initial 
condition to the stationary state, or due to a finite-speed, 
or nonadiabatic, driv ing . In any cas^ S(^) is known as 
the nonadiabatic |36l438l| or excess [30|, l39| entropy pro¬ 
duction, in contrast with the entropy production needed 
to maintain the stationary state, which is often referred 
to as adiabatic or house-keeping entropy production . 

The fluctuation theorem stated in (EHl) exploits the 
dynamical symmetries of the process through the dual 
map and the nonequilibrium potential, in the same spirit 
as the detailed fluctuation theorem for processes con¬ 
necting nonequilibrium states developed by Esposito and 
Van den Broeck Finally, the integral theorem 

(HZl) is the quantum version of the Hatano-Sasa theo¬ 
rem [ 3 ^, extending the Jarzynski equality to nonequilib¬ 
rium states. The corresponding second-law-like inequal¬ 
ity da extends to arbitrary boundary conditions the 
quantum Hatano-Sasa inequality for concatenated CPTP 
maps proposed by Sagawa [^. 

IV. APPLICATIONS 

Despite their simplicity, the above fluctuation theo¬ 
rems include as special cases many of the known quan¬ 
tum fluctuation relations. In the section, we explain how 
these relations come about in our formalism. We first 
discuss the boundary term a(n,m) and then apply the 
general theorem to different dynamics. Here we specify 

= 0, the anti-unitary time-reversal operator. 


A. Boundary terms 

There are two common choices for boundary terms: ij 
setting the initial state of the dual equal to the final state 
of the forward process pf = ApiA ^, pi being an arbitrary 
state; and ii) setting the initial state pi of the forward 
process and the initial state pf of the dual process as 
equilibrium states. Notice that by selecting the initial 
states of the forward and dual processes we are also fix¬ 
ing the basis in which the quantum measurements are 
performed at the beginning and end of the processes. 

In the first case, the boundary term 

<y{n, m) = — Inpf (m) -|- Inpi(n) = Sf (m) — Si(n) (22) 

is the increase of the stochastic or trajectory entropy fl^ . 
nil. Hi], 113, whose average over forward trajectories yields 
the increase of von Neumann entropy. 

This choice is relevant from a theoretical point of view, 
but the resulting dual process is hard to implement in 
general, except when the system is small enough to be 


prepared in an arbitrary state (say, a few qubits or a 
harmonic oscillator). 

The second choice, equilibrium initial states for the 
forward and dual dynamics, is more interesting from an 
operational point of view, since the dual dynamics can 
be easily implemented in the laboratory by equilibrat¬ 
ing the system with a thermal reservoir and reversing 
the protocol that drives the Hamiltonian d, [l^, . Let 

us suppose that, before applying any quantum map, the 
system Hamiltonian is initially fixed iLi, whereas after 
the Hamiltonian is Hf. We further take the initial state 
of the forward process to be equilibrium at inverse tem¬ 
perature /3, that is. Pi = where Fi is the cor¬ 

responding free energy. Similarly, we initialize the dual 
process in the final equilibrium at the same temperature, 
pf = Then, 

a{n,m) = p(El-El-Ff + Fi) = /3(AE;„,^-AE) (23) 

where the are the eigenvalues of the initial and 

final Hamiltonians, respectively. 


B. Unital work relations 

As a first example, we take our quantum map to be 
unital (or bistochastic 0), that is, the identity is an 
invariant state, f(Il) = 1 (although the identity may 
not be the only one). Any unitary evolution U is uni¬ 
tal, = 1 , and its dual map is the time-reversal 

U = 0[/^0^. Another example of a unital map is the 
projective measurement of an observable but, more gen¬ 
erally, any minimally disturbing measurement is unital 
(0 . For these maps, the Kraus operators are self-adjoint 
= Mk, leading to dual operators Mk = 0Mfe0l. Fi¬ 
nally, pure decoherence is also implemented with unital 
maps that remove all the off-diagonal elements in a spec¬ 
ified basis. For all such unital maps or concatenation of 
such maps, A$,r(fc) = 0 for all k, and the fluctuation 
theorem only consists of the boundary term. 

Let us now consider a concatenation of unital maps 
as describing a physical process. An important exam¬ 
ple is a process consisting of several unitary transfor¬ 
mations induced by driven time-dependent Hamiltoni¬ 
ans, punctuated by a number of measurements and/or 
pure decoherence processes. In each map, energy can 
be transferred to the system. We call the energy input 
into the system due to the driving Wdrive, driving work, 
and Wmeas the energy input due to the measurements 
and/or decoherence processes. Whereas the driving work 
Wdrive has & dear interpretation as the energy supplied 
by driving, the origin of the energy input due to mea¬ 
surement is still obscure. This energy transfer occurs, 
for instance, in a projective measurement of an observ¬ 
able that does not commute with the Hamiltonian. In 
any case, AEn,m = Wdrive + Wmeas and, if we choose 
equilibrium initial states the boundary term a is given 
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by (l23l) and 

S(7) = /3(H^drive + " AF) = /SW^diss- (24) 

The fluctuation theorem (EOl), therefore, reproduces the 
work fluctuation theorems for unital processes derived in 
[Hill, [Hill (see also [i5|,|4^). Notice that, if we allow 
the system to relax to equilibrium after the maps have 
been applied, then 11 ( 7 ) equals the entropy production 
along the whole process. We stress that this result is valid 
for any concatenation of unital maps. On the other hand, 
if we choose the initial state of the dual process as the 
final state of the original process, E = —\npf{m)/pi{n) 
is just the change in stochastic entropy. When averaged, 
the entropy production E becomes the change in the von 
Neumann entropy of the system 

AS-sy, = (E) > 0, (25) 

whose positivity follows from (flSl) . This provides another 
proof of the well-known property that unital maps can 
only increase the von Neumann entropy 0. 

C. Thermalization and heat 

Another interesting example is a generic thermaliza¬ 
tion map [T^ at inverse temperature /? = l/{kT) (or 
Gibbs-Preserving map (T7I|1. that is, a map whose invari¬ 
ant state is the equilibrium density matrix tt = 
where H = Ej \ej) (e^j is the Hamiltonian of the sys¬ 
tem and F its free energy at temperature T. Thus, 
the nonequilibrium potential is related to the energy as 
<i)^(j) = — ln7r(j) = /3(F — Ej). To verify our fluctu¬ 
ation theorem, each Kraus operator Mk must promote 
transitions between energy eigenstates involving a given 
change of energy AE{k), that is, Mk = Y.ji \^j) (siL 
where the sum runs over pairs of energy eigenstates with 
the same energy difference AE{k) = Ej — Ei. Now, since 
the energy is supplied by a thermal reservoir, we can 
identify these energy exchanges as heat flowing to the 
reservoir, Q{k) = —AE{k). The dual Kraus operators 
Mk oc mI = J2ji \^i) (Cj I (for a time-reversal invari¬ 
ant H) induce the reverse transitions accompanied by 
the reverse flow of heat Q{k) = —Q{k), and thus can be 
identified with a Kraus operator in the original map. 

We now can consider a thermodynamic process formed 
by a concatenation of thermalization steps by N distinct 
thermal reservoirs with inverse temperatures {I3i}^i in¬ 
terspersed by unital transformations (unitary drivings, 
measurements or decoherence). For this setup, if we 
choose the initial state of the dual process as the final 
state of the original process, we arrive at 

N 

E(7) = Sf(m) - Si(n)-b^AQi(7), (26) 

i=l 

with Qiij) the total heat flow into the i-th reservoir. 
In this case, we get a fluctuation theorem for the total 
irreversible entropy production in the process. 


On the other hand, the equilibrium boundary terms 
are interesting when restricted to one thermal reservoir, 
leading to 

E(7) = /3(AA„,^ - AA + Q(7)) = - AF), (27) 

employing the energy balance AEn^m = bF(7) — Qil)- 
Again, E equals the entropy production along the whole 
process consisting of the map concatenation followed by 
a thermal relaxation. The detailed and integral fluctua¬ 
tion theorems following from the identification (l27l) are 
respectively the quantum Tasaki-Crooks and Jarzynski 
fluctuation theorems for thermal maps punctuated by 
unital maps 

D. Lindblad master equations 

Another nice illustration of our results are the Lind¬ 
blad master equations that model the Markovian dy¬ 
namic evolution of open quantum systems Bm. For a 
quantum system with Hamiltonian iJ, a Lindblad master 
equation is specified by a collection of positive Lindblad 
operators as 

dtpt = -i[H, pt] + 'E[Lk]pt = Cpu (28) 

k 

where the dissipator V is defined as T’[L]p = LpL^^ — 
^ Lp + pL^Ly To make contact with our fluctuation 
theorem, we begin by observing that the solution to (|2^ 
can be obtained by concatenating a sequence of maps 
together that evolve the system forward in small time 
steps dt: 

K 

^{Pt) = (1 + Adt)pt = MqPiMq + MkPtM^., (29) 

fc=i 

with Kraus operators 

Mo = (30) 

k 

Mk = LkV^, 1 < k < K. (31) 

This map has at least one invariant state tt [ 3 ^, obeying 
Lit = 0. 

To satisfy our fluctuation theorem, the Kraus opera¬ 
tors {Mk} must be of the form (fT^ and verify the gen¬ 
eralized detailed balance relations da. Enforcing these 
conditions on {Mk}k>i immediately leads to a restric¬ 
tion on the Lindblad operators similar to (na- Namely, 
each Lindblad operator must induces jumps between 
invariant-state eigenstates, Lk = Kj) where 

my = 0 for all i, j such that — 4>^(z) 7 ^ A^T^{k). In 

this case, the generalized detailed balance relation m 
holds: 

Lk = eLley k>l. 


(32) 
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As for the Kraus operators, if the Lindblad operator 
induces jumps where the nonequilibrium potential 
change equals A$ 7 r(fc), then they obey commutation re¬ 
lations similar to (TT^ : 

[ife,ln7r] = 

(33) 

and [Lj)Lfc,ln7r] = [L\Lk,TT\ = 0. 

Let us verify now whether Mq also satisfies our condi¬ 
tions. The dual operator ([81) reads: 


Mn = 07r2 


(34) 


Since [Lj)Lfc,7r] = 0, for our generalized detailed balance 
condition to hold, that is Mq oc 0M^0f, we must assume 
that [iL, tt] = 0, forcing the invariant state to be diagonal 
in the energy eigenbasis. An immediate consequence of 
this observation is that in fact the A^.^- must correspond 
to jumps in the energy. With this additional assumption, 
we have 


Mo = ©Ji- 


= 0Mf|0t 


(35) 


Thus, Mq satisfies our generalized detailed balance re¬ 
lations with A<h 7 r( 0 ) = 0, as one would expect for a 
Kraus operator that does not induce transitions. The 
restrictions on the Lindblad operators outlined here as 
assumptions can in general be proved as consequences of 
the requirement that the dual map also be CPTP 

Consider now the following process. We run the Lind- 
bladian evolution for an interval of time [ 0 ,r], and mea¬ 
sure some observables at time t = 0 and t = t. In this 
scenario, a trajectory 7 = (n, fci, /c 2 ,..., m) is given 
by the initial and final measurement outcomes, n and 
m respectively, and a set of jumps ki occurring at times 
ti. Notice that the stochastic trajectory, as defined in 
the previous sections, should contain a big number of 
instances k^. = 0 , i.e., corresponding to operation Mq, 
between jumps. However, these operations do not con¬ 
tribute to 51 ( 7 ) and we can omit them from the discus¬ 
sion. In this case. 


S( 7 ) = cr(n, m) + ^A<i>.(fc,), (36) 

i 

With the entropic boundary conditions (l2^ . we arrive 
at the quantum generalization of the Hatano-Sasa theo¬ 
rem [ 3 ^ for the nonadiabatic entropy production of Lind¬ 
blad master equations, as developed in [l^. Its average 
over trajectories, corresponds to the integrated expres¬ 
sion first introduced bySpohn for arbitrary quantum dy¬ 
namical semigroups |49| . then extended by Yukawa to 
driven quantum Markov processes Q. The equivalence 
between our trajectory picture and the average thermo¬ 
dynamics behavior has been discussed in Q. 


So far we have been treating the dissipation in the 
Lindblad master equation as a whole. When the dissipa¬ 
tion can be interpreted as coming from M distinct ther¬ 
modynamic reservoirs (or Markovian noise processes), we 
can employ our formula for the entropy production of 
concatenated maps dini) to arrive at a complementary 
formulation of the thermodynamics. The effect of each 
of the M reservoirs is captured in the dynamics by a sep¬ 
arate collection of Lindblad operators where 

a = 1 ,..., M labels the reservoir: 

dtPt = -i[H,Pt] + EE V[Lk,c.]pt. (37) 

a k 

Similar to before (1291) . we can implement the evolution 
of this equation over a small time interval dt by a map, 
except now it is formed by a concatenation of inter¬ 
mediary maps, £{pt) = £aM ''' ^ai^oipt), each arising 
from the different dynamical influences. The first map 
£o{Pt) = Pt — 'i[H^pt]dt captures the unitary part of the 
dynamics with a single Kraus operator = 1 — iHdt; 
the subsequent maps describe the dissipative reservoirs, 
whose Kraus operators are 

Mo,a = 1-(38) 

k 

Mk,a = I < fc < Ka- (39) 

Notice that the exact sequence of maps £a is immate¬ 
rial as they all commute to first order in dt. Crucially, 
each reservoir is assumed to have its own invariant state, 
^ 0 ( 7 ’"^“^) = 7r(“) (or equivalently = 0). 

For example, a thermal reservoir at inverse temperature 
/3(“) would have the equilibrium Boltzmann density ma¬ 
trix 7r(“) = as its invariant state. The corre¬ 

sponding Lindblad operators must then induce jumps in 
that state, Lk,a = ^ satisfy our 

generalized detailed balance relation dig. As a result, 
the {Mo,a}^o immediately satisfy the generalized de¬ 
tailed balance relations with A<I>= 0, which remark¬ 
ably does not require the invariant state to commute with 
the Hamiltonian. 

Now, a trajectory for this setup corresponds to a list 
7 = (n,k^^\k^‘^\ ... ,k^^^\m)^ given by the initial 
and hnal measurement outcomes, n and m, and a set 
of jumps occurring at times U in the ai reservoir. 
Notice that only one jump in one of the M reservoirs can 
happen in any given dt, since the probability to observe 
two jumps is negligible. The result from (1^ is then 

S( 7 ) = cr(n, m) 4- ^ (A:}“‘^). (40) 

i 

This point of view allows us to treat multiple reservoirs at 
once, such as an engine operating between a hot and cold 
thermal reservoirs, each represented by a different set of 
Lindblad operators [5^. Using the entropic boundary 






conditions (I22p . the resulting average entropy produc¬ 
tion has l ong been known from the works of Spohn and 
Lebowitz [51| and Alicki et al. [H^- 

Remarkably, condition (155)) is fulfilled by almost all 
known examples of driven Lindblad equations for systems 
weakly coupled to reservoirs. If the Hamiltonian H of 
the system is constant, the weak coupling limit results 
in a Lindblad equation where the operators are 

labelled by the Bohr frequencies ui which are transition 
frequencies between the levels of the Hamiltonian H, i.e, 
they are of the form lo = uii — ujj, for some pair of levels 
i, j with energies = huJi and ej = tvjjj, respectively [ 1 , 
l34| . These are ladder operators that lower and raise the 
energy levels of H, obeying the commutation relations: 

[L^,iL]=ccL^ ; [Li,H] = -ujLl (41) 

Their commutator with the logarithm of the stationary 
density matrix can be written as: 

( 7 ri|[L<^,ln 7 r]| 7 rj) = ( 7 ri|L^| 7 rj) In (42) 

For dMl) to be satisfied it is sufficient that the ratio 
Tr{i)/Tr{j) = is a function of the energy differ¬ 

ence Acij = Cj — Ci- In that case 

[L^,ln7r] = f(nuj)L^ (43) 

and A(/)^(a;) = /(feu). In the case of a single thermal 
reservoir /(e) = fie, and A(f>T^{u}) is entropy production 
in the reservoir (heat flow divided by temperature) asso¬ 
ciated to a transition of frequency w. Furthermore, the 
Lindblad operators will come in pairs such 

that To; = L-uj oc LJ,, and every jump can be undone. 
As a result, the dual process is equivalent to the original 
process. This approach was developed for work fluctua¬ 
tions theorems in and heat fluctuations in [s^ . 

The preceding arguments can naturally be extended to 
a time-dependent Hamiltonian Hft) and time-dependent 
Lindblad operators Lk{t), yielding an instantaneous sta¬ 
tionary state 7 r(f) (or states 7 r(“)(t)) [13, HI]. This is 
the case when the Hamiltonian H{t) = H{Xt) is driven 
through the slow change of a collection of external pa¬ 
rameter At, the Lindblad operators become parameter¬ 
ized by the external parameters Lfc(At), and our gener¬ 
alized detailed balance relation will hold at every time 
d, 113 . They even continue to hold for nearly adiabatic 
driving . For fast periodic driving, Floquet theory can 
be used to derive a Lindblad master equation [sl]. This 
theory picks out as a preferred eigenbasis a collection of 
time-periodic states, or Floquet states, each with a corre¬ 
sponding quasi-energy or Floquet energy. The collection 
of Lindblad jump operators {Lfc} then induce transitions 
between Floquet eigenstates of the periodic Hamiltonian 
leading again to the generalized detailed balance relation 
(1551) with A$,n.(fc) the change in Floquet eigenvalues in 
the fc-th jump, which often corres pon ds to the heat ex¬ 
hausted into the environment |33l. l57l|. Finally, our pre¬ 
dictions can be used to recover the fluctuation theorems 
derived for driven Markov dynamics presented in [l5l| . 


It is remarkable that our fluctuation theorem can yield 
different results for S, depending on the resolution of the 
stochastic trajectory. For instance, in the case of the sys¬ 
tem in contact with several thermal reservoirs, E is given 
by (l40l) if the trajectory keeps track of the jumps induced 
by each reservoir separately. On the other hand, if the 
trajectory only gives information about the jumps of the 
system in the basis where the stationary density matrix 
of the entire Lindblad equation is diagonal, we have dMl). 
Consequently, for the same map one can have both (1551) 
and (HO)) . The distinction is the same as the difference be¬ 
tween the fluctuation theorem for the entropy production 
(l40l) and the nonadiabatic entropy production (156)) . 


V. CONCLUSION 

We have presented a general fluctuation theorem for a 
large class of completely-positive, trace-preserving quan¬ 
tum maps that verify the generalized detailed balance 
condition in on. From these relations many of the 
known quantum fluctuation theorems follow naturally. 
Included in this family are classical fluctuation theorems 
for arbitrary stochastic maps, as such maps are special 
cases of CPTP quantum maps where the dynamics re¬ 
main diagonal in a particular basis. The theorem exploits 
the dynamical symmetries of a process and its dual and 
can be interpreted as a quantum version of the Hatano- 
Sasa theorem [s^. When specialized to maps induced 
by thermodynamic reservoirs, our results reproduce the 
quantum fluctuation theorem for entropy production. 

We have extended the notion of the dual process, intro¬ 
duced by Crooks [l3| and clarified its relationship with 
the time-reversal process used by Campisi et al. and 
Watanabe et al. to derive fluctuation theorems for uni- 
tary evolution punctuated by projective measurements 
[3, [l3 and with the classical dual process used by Espos¬ 
ito and Van den Broeck to split the entropy production 
into an adiabatic and nonadiabatic contribution [s^. 

For nonunital maps our work should be contrasted with 
the integral fluctuation theorems presented in [ 3 , , 

which in our notation reads 

(e”'") = 7, (44) 

for a process dependent correction factor 7 . We have 
seen that by including the nonequilibrium potential, no 
correction term is necessary The resulting E, can then be 
given a clear interpretation as an entropy production in 
most setups of physical interest. 

Our results also show the peculiarity of unital maps 
regarding entropy exchange, as already pointed out in 
[ 13 . [l3 . [l7l - [l9j . The nonequilibrium potential associated 
to those maps is constant and therefore it does not ap¬ 
pear in the fluctuation theorem. The entropy production 
E in this case is only given by the boundary terms, sug¬ 
gesting that unital maps can be induced without any en¬ 
tropy exchange between the system and its surrounding. 
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Thermalization at infinite temperature is an obvious ex¬ 
ample, but decoherence or, equivalently, projective mea¬ 
surements, are relevant examples of unital maps. In all 
these cases, energy exchange between the system and its 
surroundings can occur, but this energy exchange does 
not imply any entropy change in the environment. 
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